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ABSTRACT 

We present high-resolution infrared spectroscopy of the Herbig Ae star HD 141569 A in the CO v = 
2 — 1 transition. With the angular resolution attained by the adaptive optics system, the gas disk 
around HD 141569 A is spatially resolved down to its inner-rim truncation. The size of the inner 
clearing is 11±2 AU in radius, close to the gravitational radius of the star. The rough coincidence to the 
gravitational radius indicates that the viscous accretion working together with the photoevaporation 
by the stellar radiation has cleared the inner part of the disk. 

Subject headings: circumstellar matter — planetary systems: formation — planetary systems: proto- 
planetary disks — stars: formation — stars: individual (HD 141569 A) — infrared: 
stars 



1. INTRODUCTION 

Direct imaging of circumstellar disks is one of the most 
rapidly developing fields of observational astronomy (e.g. 
Grady et al. 2005). It faces, however, two barriers in the 
optical/infrared: it is mostly sensitive to the dust grains 
either in scattered light or thermal emission, although 
dust carries only a fraction of the mass with all the rest 
in the gas phase mostly in molecular form. Moreover, it 
probes only outer regions of a disk, being obstructed by 
glaring light from the central star. However, it is this 
invisible part of the disk, gas inside 30 AU of the star, 
that has direct consequence to the gas giant planet for- 
mation. The gas dissipation at the inner disk practically 
limits its time scale, and make a critical test to the cur- 
rent planet formation theories that predict distinct time 
scales to build up Jupiter-like planets (~10^ yr. Pollack 
et al. 1996; -10^ yr. Boss 2002). 

Molecular spectroscopy in the near-infrared has a par- 
ticular advantage in this area. The energy gap between 
the ground state and the vibrationally excited states is on 
the order of 10'^ K, which automatically guarantees that 
the gas in emission is in the innermost part of the disks. 
The molecular spectroscopy also ensures that the emis- 
sion does not arise as close as in the stellar photosphere 
where the temperature is too high to keep the molecules 
from dissociation. The CO molecule is the versatile spec- 
troscopic tracer of the gas disks use d by many authors 
both at the first overtone at 2.3 um lIGeballe fc PerssonI 
19871 ICarriligsatlChandler. Carlstrom. fc Scovillel | 199.- ' 
"" 596';' 



Naiita et a l."1996:'Biscava et al.'1997*.'Kraus et al j l2f 
Bik fc Thi 2004; Blum ct al. 2004; Thi et al. 2009), aj^d 
the fundamental band at 4.7 A tm (Brittain et al. 20031 
[Naiita. Carr fc Mathieul 120031: iBlake fc Boogert 20041 
H^ttig ct al.l l2004HCarmorna et alJl2005() . 
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HD 141569 A is a Herbig Ae star in transition to a de- 
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bris disk object with an optically thin disk with relatively 
small infrared excess (Lir,/L* = 8 x 10"'^ by Sylvester et 
al. 1996; cf. Lir/L* = 2.4 x lO'^ of (3 Pic by Heinrich- 
sen et al. 1999). The dust disk around HD 141569 A 
has been directly imaged with coronagraphy by HST 
in the scattered light at the visible and the near- 
infrared wavelengths (Augereau. Lagrange. & Mouille. 
n"999: W einberger ct al. 1999 ; Mouillct ct al. 20q1 
IClampin et aLl2003j) . as well as a ground-based tele scope 
with adaptive optics system ijBoccaletti et al.l2003j) . The 
position angle of the dis k is 356°±5° east frorn north 
in projection on the sky ({Weinberger et alJll999() . The 
disk ellipse in the thermal emission is well aligned to 
the dust scattered hght (355° ±19° by Fisher et al. 
2000; 354°±4° by Marsh et al. 2002). The resid- 
ual gas in the outer disk has been detected in the ra- 
dio spectroscopy at CO J = 1 — and J = 2 — 1 
(Mh2 — 300 Mq, Zuckerman, Forveille, & Kastner 
1995; Dent, Greaves, & Coulson 2005). The line emis- 
sion CO J = 2 — 1 is recently spatially resolved by 
iDutrev. Lecavelier Des Etangs. fc Augereaul (|2004) with 
red- and blue-shifted components split in north and south 
of the star with a velocity interval indicative of gas or- 
bital motion. 

HD 141569 A is exceptional among Herbig Ac/Be stars 
for its CO Tj = 2 — 1 emission. The infrared transition 
V = 2 — 1 manifests the molecul es are first electroni- 
cally excited by the UV irradiation. iBrittain et al.1 l|2003D 
argues that the absence of high velocity wings in CO 
u = 2 — 1 indicates that the line emission comes from the 
inner-rim of the disk receded from the star by 17 AU, 
where the disk wall is directly illuminated. The pres- 
ence of inner clearing of the size of 10-20 AU is consis- 
tent with the lack o f near-infrare d excess in the infrared 
energy distribution ijSvlvester et al 19961. and the mid- 
infrared imaging implyi ng disk withi n 30 AU is already 
cleared to some extent l|Marsh et al.ir2002^) . The size of 
the inner hole is well within the reach of spatially re- 
solved observations with adaptive optics systems at 8-m 
class telescopes at the distance of the star [d =99-108 pc, 
van den Ancker et al. 1998; Merin et al. 2004; hereafter 
we take d =108 pc from Merin et al). The goal of this 
paper is to present such directly resolved observations of 
this inner clearing in CO w = 2 — 1 transitions. 
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2. OBSERVATION AND DATA REDUCTION 

The spectroscopic observation was carried out on 
25 May 2005 UT at the S ubaru Telescope with 
the facihty spectrogra ph IRCS ijTokunaga et al.l Il998t 
iKobavashi et alJ 1200(H) • The curvature sensing adap- 
tive optics system was used to feed nearly diffraction- 
limited images to the spectrograph (Gacsslcr ot al. 2002; 
iTakami et al. 2004). HD 141569 A (i?=7.02; Meri'n et al. 
2004) itself was used as the wavefront reference source 
at the visible wavelength. The grating angles were set 
so that the two echelle bands at 4.597-4.718 /im and 
5.014-5.147 fj,m were covered in the Ikxlk InSb array. 
The first short-wavelength band includes the i?-branch 
of CO t; = 2 - 1 at R{0) to i?(15). The resolving power 
of R =20,000 (Aw =15 km s^^) was attained using a 
narrow slit of 0'.'15x7'.'5. The plate scale along the slit 
is 60 mas per pixel. The slit was oriented to north to 
south, parallel to the major axis of the disk on the sky. 
The data were recorded by switching the field of view by 
3" in every 5 minutes along the slit with the tip-tilt mir- 
ror inside the adaptive optics system in order to subtract 
the sky emission. The total integration time is 20 min- 
utes. The spectroscopic flat field was obtained at the end 
of the night from a halogen lamp exposure. The seeing 
was modest, 0'.'9 at R, during the observation. 

Two-dimensional spectrograms were reduced in the 
standard manner with the IRAF^ image reduction pack- 
age, involving sky-subtraction, flat-fielding, and interpo- 
lation of outlier pixels. The wavelength calibration was 
carried out by maximizing the cross-correlation between 
the observed spectra and the atmospheric transmission 
curve calculated by ATRAN (Lord 1992). The tih of the 
dispersion axis with respect to the array column, and the 
geometrical distortion, i.e., the curvature of the slit image 
on the detector array, was corrected so that the disper- 
sion and the telluric emission lines show up straight to 
the detector vertical and horizontal over whole spectral 
range. A part of the spectrogram after the correction of 
spectral and spatial distortions is shown in Figure ^ 

The CO line emission at t; = 2 — 1 within our spec- 
tral coverage is all found spatially extended up to 50 AU 
beyond the photospheric emission of the star. In or- 
der to isolate the molecular lines from the stellar photo- 
spheric emission, one-dimensional point spread function 
(PSF) was taken from the nearby continuum, and was 
subtracted from the line emission. Note that the spec- 
tral PSF, although spatially one dimensional, provides 
almost an ideal PSF: it is recorded by the same instru- 
ment at the same time at nearly identical wavelengths 
to the line emission. The size of the PSF measured at 
the continuum wavelength was 0'.'24 in full width at half 
maximum. Caution was taken at the scaling of the PSF, 
since over-subtraction of the PSF introduces an artificial 
depression at the center, which is possibly mistaken for 
a disk inner hole. First, one dimensional spectrum was 
extracted inside a small aperture matched to the PSF to 
sample the stellar continuum flux exclusively. The one 
dimensional spectrum was smoothly interpolated from 
the both sides of a line to measure the continuum con- 

® IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 




Fig. 1. — Left: a cut-out of the spectrogram of HD 141569 A 
at 4.6 fim obtained with IRCS at the Subaru Telescope on the 
night of 25th May 2005. The slit was oriented from north to south 
along the major axis of the disk projected on the sky, with north 
being the left hand side. The spatial and spectral distortion was 
corrected using emission lines in the telluric atmosphere. The emis- 
sion lines of CO v = 2 — 1 within our spectral coverage are all found 
spatially extended up to 50 AU at both sides of the star. Middle 
column: the blow-up of the emission lines from R{5) to R{9). Right 
column: same as in the middle panels, but after one-dimensional 
point spread functions have been subtracted. The centrer of the 
PSF, therefore, the position of the star is marked by dotted verti- 
cals. Disk rotation is clearly visible, with the northern disk reced- 
ing from us. The rota tion is consistent wit h the spiral arms seen 
in the dust scattering tClamp in et al.|l2003). and the earlier radio 
spectroscopy (Dutrey et al. 2004; Figure 3 presented in their paper 
accidentally shows the northern component of the CO emission is 
approaching toward us. The confusion has been cleared by private 
communication with A. Dutrey and J.-C. Augereau). 



tribution at the line wavelength. The PSF was scaled to 
the height of the continuum flux, and subtracted from 
the line emission accordingly (Figure^. 

We briefly discuss the influence of the telluric absorp- 
tion below. The interference of CO emission lines over 
the telluric absorpti on is illustrated in the spectra in 
iBrittain et al.l l)2003j) obtained with similar spectral res- 
olution with the present observation. The energy loss in 
absorption in the telluric atmosphere has to be compen- 
sated to restore the absolute flux in the line emission, and 
pottentially the source of large systematic uncertainty in 
the flux measurement. The radial profile is, on the other 
hand, less subject to the atmospheric effect, as long as 
the emission line stays at constant wavelength along the 
radius. However, because of the finite displacement in 
the line velocity by the gas orbital motion, the emis- 
sion lines still have to be at reasonably flat part of the 
transmission curve. Fortunately the interference of the 
telluric absorption is less severe at CO v = 2 — 1 than 
COu = l — Oin which the telluric atmosphere has corre- 
spond ing absorption lines, as is also seen in Brittain et alJ 
l|2003() . Since apparent gas velocity at 5 AU and 20 AU 
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is 15 km s^^ and 7 km s^^ respectively, if it is in the 
Keplerian rotation, a flat part of Av « 8 km s^^ would 
suffice for the radial profile to be intact. Those lines 
longword of 4.68 fim that seriously overlap with the tel- 
luric absorption lines with the clearance less than this 
interval were discarded, and the lines from i?(5) to -R(9) 
were only used for the further analysis at reconstruction 
of the radial profile and the line image shown in Figure |2 
and El 

3. RESULT 

Each line image was summed up together along the 
dispersion axis by ±19 km to reconstruct a one- 
dimensional disk radial profile. The central depression 
is clearly resolved in the combined profile in Figure [3 
Since the individual line image is considered as an inde- 
pendent imaging, the error bars were given by the stan- 
dard deviation of the radial profiles at different transi- 
tions. A power-law radial profile with an inner cut-off 
at ri {f{r) oc at r > r^, /(r) = at r < r^) was 
fit to the observation after convolved with the observed 
PSF. The best-fit radius of the central clearing is found 
at Ti =11±2 AU with the power-law index /3 =2.6±0.6. 
The combined disk image with improved pixel sampling 
is shown in Figure |31 Disk rotation is clearly detected 
as a velocity offset between the southern and northern 
disks, approaching and receding from us, respectively. 
HD 141569 A has its eastern disk near to us known 
from the asymmetric azimuthal illumination of the disk, 
which is attributed to the forward scattering of the grains 
(^cinbcrgcr et al. 1999). The disk rotation is therefore 
clockwise in projection on the sky, which is consistent 
with the winding of the sp iral arms at 200-400 AU of the 
star l|Clamnin et alJl20?)3) . A simple geometrical model 
is calculated for an infinitely thin disk vignetted by the 
slit field-of-view with the gas in the disk in the Keple- 
rian rotation until the inner truncation (M, = 2.0Mq, 
Marin et al. 2004; inclination angle (j) = 51° from face-on, 
Weinberger et al. 1999). The velocity dispersion along 
the line of sights well agrees with the observed disk image 
in the line emission, and the lack of high velocity com- 
ponent close to the central star manifests the innermost 
part of the disk is indeed evacuated (Figure ISJ. The 
line emission from the northern disk in Figure 3 looks 
slightly brighter than that from the southern disk. The 
asymmetry might be also recognized in each line pro- 
file shown in Figure 1 before it is combined together, 
however, more observations are necessary to confirm the 
asymmetry, and discuss further on the possible cause of 
it. 

4. DISCUSSION 

The size of the central cavity gives some insights 
into its formation mechanism. The inner clearing is 
apparently larger than the magnetospheric truncation 
(r ~ i?* ~ 10~^ AU), or the dust sublimation radius 
(r ~ 0.1 AU), therefore those mechanisms are not likely 
responsible. The cavity size is comparable to the gravi- 
tational radius of the star where the sound speed in the 
ionized medium is equal to the escape velocity of the sys- 
tem (rg = GAh/cf; Shu et al. 1993, HoUenbach et al. 
1994; «18 AU for Af, = 2.0 Mq). The gravitational 
radius defines the innermost radius of a disk where the 
gas is driven away from the system once ionized by the 
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Fig. 2. — The radial profile of CO line emission at the central 
part of the disk. A power-law radial profile (dotted line) was fit to 
the observation after convolved with the PSF (solid line) to find 
the best-fit cut-off radius at = If it 2 AU. The error bars (f cr) 
are given by the standard deviation of the radial profiles extracted 
from CO ?; = 2 - I R{5) to _R(9). 



stellar radiation. The rough agreement of the inner edge 
of the disk around HD 141569 A with the gravitational 
radius of the star indicates the cavity is being cleared by 
the photoevaporation working in combination with the 
viscous accretion. 

The dissipation of an inner circumstellar disk is 
controlled by the two processes: the mass loss by 
photoevaporation and the viscous accretion on to the 
star. The viscous accretion is the primary drive to 
remove a disk, transporting bulk of the material in the 
outer disk inward eventually onto the star. On the 
other hand, the photoevaporation is a slow process, 
steadily removing a disk outside the gravitational radius 
as long as the extreme ultraviolet irradiation (EUV; 
hv > 13.6 eV) maintains. However, the persistency or 
EUV radiation, or, the origin of the ionizing photon is 
unclear, whether it is a consequence o f active accretion 
ifMuzerolle. Calvet. fc HartmannI 1200 H iMuzcroUe et alJ 
12004^, or of the stellar activity (e.g. Deleuil et al. 
2005). In the former case, the photoevaporation is 
expected to drop off as the viscous accretion slows 
down, whi ch causes problem to disperse a disk within 
10'' y r s fMatsuv ama. Johnstone, fc HartmannI '2003; 
iRudenI 12004 ) . Although the photoionizing radiation 
is likely chromospheric, and holds u ntil the end of 
the d isk evolution in low-mass stars l|Alexander et alJ 
120051): it may not directly apply to Herbig Ae/Be stars 
iCxradv et all200l . The mass loss via photoevaporation 
proceeds most effectively immediate outside of Vg, as 
the disk wind decreases with the radius as S(r) oc r^^/^ 
l|Hollenbach et all Il994|) . During the early phase of 
the disk evolution, the photoevaporation is relatively 
unimportant, as the viscous accretion quickly fills in 
the mass loss at the gravitational radius. However, as 
the disk declines, the viscous accretion slows down, and 
is eventually balanced by the photoevaporation. The 
net inflow across the gravitational radius is quenched, 
as the photoevaporation virtually reverses it. Starved 
for the further supply, the material inside rg is quickly 
drained onto the star on the time scale of viscous 
evolution. In consequence, a disk dissipates on two dis- 
tinct time scales (Hollenbach. Yorke. & JohnstoM 

Clarke. Gcndrin . fc Sotomavoii 
lArmitaee. Clarke, fc Fallal ^ 12005 " 
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Johnston et al] I200I iTakeuchi. Clarke, fc LinI l2005t 
Alexander et al.ll2nO(iD : the whole disk lifetime may last 
as long as ~ lO"^ yrs, while the evacuation of the inner 
disk is almost an instant (At ^10^ yrs). 

The presence of the central cavity clears several issues 
how the inner disk is removed at the end of the disk 
dissipation. The photoevaporation and the viscous ac- 
cretion is likely effective to remove a inner disk, inside 
outward, without calling for any hypothetical planets to 
sweep it up. The radius of the inner clearing measures 
= 11 AU, slightly smaller than the nominal gravita- 
tional radius of Vg = 18 AU. The smaller inner radius 
is in fact preferred by the latest models that incorporate 
far ultraviolet radiation (FUV; 6 eV< hv <13.6 eV) that 
thermally heats up the gas in the disk (DuUemond et al. 
2006 for revie w). The inass los s by photoevaporation 
formu lated by iShu et al] l)1993|) and iHollenbach et al.1 
1)1994(1 is actually the photoionizati on that assumes dis- 
sipation of the gas ionized by EUV. lAdams et al.1 ((20041 ^ 
found FUV irradiation has significant contribution to 
disk wind; and it could launch as close as 0.2-0.5 Vg of the 
central star, when the pressu re gradient at the central 
region is taken into account ijLiffmanl 1200."^ iFont et all 
12004). The location of the inner-rim of HD 141569 A 
better matches the new pictures. 

HD 141569 is a triple system with low-mass compan- 
ions B and C at 9" away. The coeva l formation of the 
syste m sets a robust age of 5±3 Myr l|Weinberger et al.1 
which agrees well with 4.7 Myr given by the full 
spectral modeling with the photometric analysis of the 
primary fjMaTnTet a l. 2004). The age of HD 141569 A 
sets a challenging time scale to form gas-giant plan- 
ets inside the gravitational radius of the Herbig Ae 
star. The conventional model of gas-giant formation 
takes 1-10 Myr to form a Jupite r-like planet by ac- 
cretion of gas ont o a rocky core Pollack et al. 
lHubickvT~B odcnhe imer. fc Lissaiierl 20051) . which is un- 
comfortably close to the standard disk lifetime (~3 Myr 
for the median disk lifetime, Haisch et al. 2001; Hil- 
lenbrand 2005 for review), although many ideas have 
been proposed to accelerate the core growth p rocess 

£Qaba, Wcthcrill, & Ikoma 2003; Alib ert et. all l2005t 
lev fc Dirksea.200&:.Klahr fc BodenheimeJl2006(l . The 
age of HD 141569 A suggests that giant planets may not 
have enough time to fully deplete the protoplanetary disk 
onto themselves, but likely have to compete with disks 
in dissipation. The subtlety in the inner disk dissipation 
and the planet forming time scale may explain a large 
scatter of disk possession among the stars in a single age 
bracket (e.g. Rieke et al. 2005), and eventually the di- 
versity of the exoplanetary systems discovered to date 
(e.g. Marcy et al. 2005). 

It is tempting to speculate on the origin of the iso- 
tope anomaly in the solar system in connection with 
the present observation. The dissociation time scale of 
CO in the normal interstellar medium is on the order 
of 10^ yrs (fcco ~ 2 X 10-1° s-^ van Dishoeck 1988), 
the molecule in the disk of HD 141569 A will be quickly 
destroyed without a shielding to the dissociating irra- 
diation from the star (A < 1100 A; van Dishoeck fc 
Black 1 988). The dust extincti on could dilute the ra- 
diation ijKamn fc Bertoldil 1200(1) . however, the extinc- 
tion cannot be too large, since the emission we observed 




Fig. 3. — Top: HST/ACS image of the debris disk around 
HD 141569 A from Clampin et al. 2003 at visible wavelength 
(A =606 nm). The central 100 AU are masked by the corona- 
graph. Middle: the image of the composite line emission from 
CO V = 2 — I R{5)-R{9). Only those lines less sever ly overlapped 
with the telluric absorption lines were used. The north is left, and 
the wavelength increases to the top. Bottom: the close-up view of 
the middle panel. The gap between the northern and the southern 
disks is clearly resolved. The velocity dispersion along the line of 
sight is calculated for a disk with an inner hole of 11 AU with the 
gas inside the disk in the Keplerian rotation (for assumed central 
mass of 2.0 Mq with the inclination angle </> =51°). The veloc- 
ity contour overlaid in the shade is consistent with the observed 
line image. The rotation curve in the case of no truncation at the 
inner disk is shown in the dashed line. The lack of high velocity 
component close to the star represents that the molecular gas is 
cleared up already at the central part of the disk. There is a possi- 
ble asymmetry seen in the brightness of the northern and southern 
disks, however, more observations are necessary to confirm it. 



here is the transition from the vibrationally excited state 
(v = 2 — 1), which needs UV photons to get pumped 
up (A <1600 A; Krotkov et al. 1980). The margin is 
not too wide. The extinction dilemma could be evaded 
by the line shielding either by molecula r hydrogen, or by 
CO itself (|van Dichoeck fc BlactdlTgM ). In contrast to 
the dust extinction that dilutes UV radiation as a con- 
tinuum, line shielding blocks the dissociative irradiation 
selectively. Self-shielding is only effective where the dust 
extinction is less than Ay = 1, but the CO is already op- 
tically thick. It is not clear if any of these conditions are 
met in HD 141569 A. The absence of near-infrared excess 
(|Svlvester et al.lll996|) only gives loose constraint, being 
consi stent with the inn er-rim as close as 0.24 AU of the 
star l|Marin et alJl200^ : although the coronagraphic and 
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mid-infrared imagin g suggests dust opacity drops at 30- 
160 A U and inward l)Weinberger et al.ll99flHMarsh et ail 
l2002j) . The column density of CO in the infrared emis- 
sion (iVco ~ 10" cm-2; Brittain et al. 2003) is too 
short to make self-shielding effective l)Lee et al.l IT? 



however, the emitting CO does not necessarily properly 
represents the column density of shielding gas along the 
line of sight to the UV source. 

With all these uncertainties, if self line-shielding is 
still at work, it may have an implication to the iso- 
tope anomaly imprinted in the early solar nebula. The 
correlated isotope fractionation in ^^O and ^^O in the 
primordial meteorites was once attributed to a super- 
nova o utburst occurr ed near the solar system at its in- 
fancy iClavton r 19931. Since then, numbers of accounts 
are proposed; i njection of nucleosynthesized materi al by 
evolved stars ((Busso. Gallino. fc Wasserburg|l200l. di- 
rect b ombardment of high energy particles l)Aleon et al.l 
120051) ■ and isot ope fractionation by self-shielding of CO 
llClavtonl l200?). When the self-shielding operates, CO 
is fractionated to C^^O as the rare isotopomers C^^O 
and C-'^^O are readily destroyed for their lower opac- 



ity, leaving more ^''O and ^'^O in atomic form. The 
isotope enrichment in the atomic gas potentially ex- 
plains the anomalous isotope ratio in the present so- 
lar sy stem (Yurimoto & Kuramoto 2 004tlLvons fc Yound 
1200 5*). Our observation found CO v = 2 — 1 emission at 
the inner-edge of the disk from the star at the distance 
where planets presumably form. We might be witnessing 
the ongoing isotope fractionation similar to that might 
happen in the early history of our solar system. 
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